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ABSTRACT 

This  paper  deals  with  the  possibilities  offered  by  the  wavelet  transform  used  in  the  analysis  of  the  transient  response  of  MIS 
circuits.  This  approach  has  certain  advantages  over  the  classical  finite  element  or  finite  difference-based  algorithms  when 
applied  to  drift-diffusion  (Roosebroeck’s)  model. 


1.  INTRODUCTION 


The  transient  response  of  MIS  circuits  is  determined  by  the  transient  response  of  the  concentration  of  carriers  in  MIS 
structures,  for  example  in  the  channel  of  the  MIS  transistors. 

At  the  most  fundamental  level  the  flow  of  the  charge  carriers  is  described  by  transport  equations  derived  from  quantum 
mechanics,  semiclassical  or  classic  models.  The  common  feature  of  these  equations  is  the  employment  of  a phase  space 
density  function  f(x,v,t)  of  the  spatial  variable  xeO,  velocity  variable  veR3,  and  the  time  variable  te(0,oo). 

A widely  used  transport  equation  is  the  Boltzmann  equation1'3: 


+ v Vx/~  V£/- Vv/  = (d,f)c 

where  U(x,t)  is  the  electrostatic  potential: 


(1) 


A U = \fdv-C 


(2) 


C is  the  doping  profile  and  the  term  id  tf)c  is  the  collision  operator. 

The  simulations  by  solving  to  the  Boltzmann  equation  are  known  to  be  very  time-consuming,  so  they  cannot  be  employed 
for  fast  circuit  simulations. 

A frequently  employed  method  to  make  the  Boltzmann  equation  more  tractable  was  found  to  be  the  ‘moment  method’ 
approach. 

By  this  method  the  hydrodynamic  model  is  derived  which  in  turn  is  reduced  to  drift-diffusion  equations,  e.g.  van 
Roosbroeck  model  which  consist  of: 

Poisson’s  equation: 


sAU  = p 


(3) 


continuity  equation  for  electrons: 


~ ~~divJn  - R(n,p) 
at  q 


(4) 


* E-mail:  kpluc@wel.wat.waw.pl 


Inti.  Conference  on  Solid  State  Crystals  2000:  Growth,  Characterization,  and  Applications 
of  Single  Crystals,  Antoni  Rogaski,  Krzysztof  Adamiec,  Pawel  Madejczyk,  Editors, 
Proceedings  of  SPIE  Vol.  441 2 (2001 ) © 2001  SPIE  ■ 0277-786X/01/$1 5.00 


323 


continuity  equation  for  holes: 


~ = — divj  + R(n,  p)  (5) 

dt  q 

where  Jn  is  the  electron  current  density: 

J„  = qMn  (UrVn  - nVU)  (6) 


Jp  is  the  hole  current  density: 

dp  = qMP(UTVP  + PVU)  (7> 

p is  the  charge  concentration: 

p=q(n-  p-C)  (8) 

R(n,p,x)  is  the  net  generation-recombination  rate,  p11=jin(n,VU,C,x)  and  |4p=Pp(p,VU,C,x)  are  the  mobilities  of  electrons  and 
holes,  respectively,  C=C(x)  is  the  doping  profile,  UT  is  the  thermal  voltage,  e(x)  is  the  electric  permittivity  and  q is  the 
elementary  charge. 

That  drift-diffusion  model  constitutes  up  to  now  the  core  of  state-of-the-art  semiconductor  device  simulation  programs.  In 
conclusion,  the  mathematical  description  of  the  transient  process  in  MIS  circuits  is  based  on  the  nonlinear  partial  differential 
equations  for  electrons  and  holes,  which  are  coupled  together  through  the  Poisson  equation.  The  crucial  element  in 
realisation  of  the  algorithm  is  effectiveness  of  this  coupling. 

2.  SOME  NOTES  ON  POISSON  EQUATION 

The  potential  distribution  is  described  using 
the  equation: 


d2u  _ Q(u,x ) 

dx 2 £ 

where : 

Q{u,x)=Q(u)+qN{x)- 
N(x)  = N0  + N](x)] 

Q,  (w)  = -q(c\  exp(-  au)  - a2  exp (au)\, 

Nq  — (<7]  — 02), 


(9) 


u is  the  potential,  Q - the  charge  density,  e - the  dielectric  susceptibility,  q - the  elementary  charge;  ai , a2  - are  dependent 
on  doping  N,  a = q/kT  [1/V],  k - is  the  Boltzmann  constant  [qV/K],  T is  the  temperature  in  Kelvin.  Introducing  the 


following  variables  in  poten-tial  units: 
rewrite  (1)  in  a somewhat  different  form 


,2  . 

' ■> 

£ 


A2 


3Zl,2 

£ 


and  the  dimensionless  variable  y=x/l  we  can 
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(10) 


< 0 when  q < 0 
> 0 when  q > 0 


Equation  (14)  is  thus  the  Helmholtz  equation,  where  the  boundary  problem  is  unambiguously  solvable  for  q>0  (where  an 
electron  charge  is  the  elementary  charge  ),  and  cannot  be  solved  unambiguously  for  q<0  (where  a pozyton  charge  is  an 
elementary  charge).  The  function  Q (u)  is  presented  in  diagram  1. 

2.1.  The  conditions  for  applying  the  finite  differences  method 

Equation  (2),  for  a network  with  constant  step  h,  takes  the  following  form: 

w,_,  - (2  + /t2n,  + uM  = h1Nj ; (11) 


By  varying  (3)  we  get: 


( 


Su,-x- 


2 + /22Q,  +h '' 


du. 


-u. 


Sui  +S  ui+]  = 0 


Fig.  1.  The  dependence  of  0(u)  l2  Is  on  u , where  q > 0 

In  order  to  ensure  that  the  solution  to  equation  (2.1)  was  stable  under  the  Dirichlet  border  condition,  that  is,  where 
du  j = SUj  — 0 , it  is  sufficient 

min|^,.I|,|<5w,.+1|)<|^/j|  8Q 

that  / \ This,  in  turn,  causes  the  dependence  K = U>0  for  each  u..  Due  to  the 

< max(j<SwM|, \Sui+]  j du 

large  gradients  of  u in  the  vicinity  y=0,  we  can  apply  the  ‘expansion  of  coordinates’  method  to  (1.1),  by  using  the 
transformation 
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(12) 


(a  + l)"-1 log(fly  + l) 

a ’ log(a  + l) 

The  coefficient  (a)  is  calculated  as  tlte  limit  of  the  series  presented  in  the  following  way: 


a. 


n+ 1 


(13) 


Where: 


0.001|(t/t)|; 

w*ff  _ v k+1f  xM^/(1+^)-Q+1. 

vvn  “ V 0 2 > 

a„ 

Equation  (2)  with  variable  “z”  can  be  written  as  follows: 


V0  is  the  begining  value. 


d2u  / Adu 
— -|og(a  + .)-  = 


(a  + 1 log(flf  + 1) 


a 


(14) 


Q(u)  + N 


r(a  + \)z  -\ 
a 


= Q[t\u,z) 


The  dependence  logio(dz/dy)  on  y is  given  in  Fig.  2. 

Equation  (6)  can  be  solved  on  a uniform  network  with  constant  step  (dz)  in  view  of  the  small  difference  in  du/dz  over  the 
entire  specimen  investigated. 

An  important  feature  of  tliis  method  is  its  accuracy,  of  evaluation  according  to  the  following  formula 


M = 


i 

\Q{x,u{x))dx 


(15) 


2.2.  The  shot  method 


Equation  (10)  is  replaced  by  the  system: 


du 

dy 


Q{u)+N{y)\ 


= UL,v 


j>=0 


= iv, ■ 


(16) 
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Fig.  2.  Dependence  logio(dz/dy)  on  y.  For  UL  = [-1,  -0.3,  1]  [V]  the  line  style:  [„ ”,  „ ”,  „•  • •”]  respectively  are 

used. 


The  conditions  for  the  stability  of  the  shot  method  are  examined.  To  this  end,  we  create  variation  equations  in  the 
surrounding  j U * (_y),  V*  (y)]  fulfilling  the  beginning  conditions 


— 5u  = Sv; 
dy 


Su{0)  = eu] 


dy  du 


Jv(o)  = ev. 


du; 


(17) 


where  €u  ,£v-  relate  to  the  modules.  Equation  (3.2)  is  written  in  the  matrix  form: 


d 

du 

"0 

l 

Su 

dy 

dv_ 

dQ 

du 

0 

u=u~ 

Sv  ’ 

Su 

V 

Sv 

y=0 

(18) 


The  necessary  condition  for  stability  of  the 
trajectory  [w  is  that  the  real  part 
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Fig.  3.  The  instability  of  the  trajectory  in  the  shot  method 


of  the  eigenvalue  of  the  matrix  are  not  positive.  Since  the  eigenvalues  can  be  expressed: 


= ±y/4 


\A,  exp(- 


au 


')+  A2  exp(aw')] 


in  the  shot  method,  any  trajectory  is  unstable,  if  its  elementary  charge  is  an  electron. 

3.  PSEUDOWAVELETS  APPROACH 


(19) 


Equation  (4)  take  the  form:  ut=Lu+Nf(u),  where  L and  N are  linear  differential  operators  and  f(u)  is  a nonlinear  function.  A 
generic  feature  of  the  solution  of  such  equations  is  that  they  behave  in  a smooth,  non-oscillatoiy  and  shock-like  manner. 
According  to  analysis  based  on  the  solution  of  the  drift-diffusion  model,  in  which  the  above-mentioned  equations  are 
coupled  through  elliptic  Poisson  equation,  the  transient  response  of  the  carriers,  when  pulsed  from  accumulation  into 
inversion,  can  be  divided  into  three  periods:  dielectric  relaxation  time,  depletion  time  and  equilibration  time.  During 
dielectric  relaxation  time,  the  surface  concentration  of  the  majority  carriers  sharply  decreases  (in  some  cases  by  about  8 
orders  in  1 ns).  During  depletion  time  it  becomes  constant  and  during  equilibration  time  it  slowly  increases  to  equilibrium 
value 

In  the  approach  presented  the  solution  u(x,t)  and  the  operators  L and  N were  projected  onto  a wavelet  base.  The  vanishing 
moments  of  the  base  functions  are  the  reason  for  the  sparsity  of  the  representation  of  both  the  solution  and  operators.  This 
sparsity  makes  it  possible  to  develop  fast,  adaptive  algorithms  for  applying  operators  to  functions  and  for  calculating  the 
functions  in  the  wavelet  base. 

These  algorithms  use  the  fact  that  wavelet  expansions  may  be  viewed  as  a localized  Fourier  analysis  with  a multiresolution 
structure  that  adapts  automatically  to  both  the  smooth  and  shock-like  behavior  of  the  solution.  In  smooth  regions,  few 
wavelet  coefficients  are  needed,  and  in  singular  regions  the  large  variation  in  the  solution  requires  more  wavelet 
coefficients.  This  approach  make  it  possible  to  combine  many  finite-differences,  (pseudo-)  spectral  and  adaptive  grid 
methods  into  a collection  of  efficient,  generic  algorithms,  named  adaptive  pseudo-wavelet  algorithms.  The  approach  was 
applied  to  transient  responses  in  some  MIS  circuits,  and  numerical  results  are  presented. 

The  approach  takes  advantage  of  the  efficient  representation  of  functions  and  operators  in  wavelet  bases,  and  updates  the 
solution  by  implementing  recently  developed  adaptive  algorithms  These  algorithms  are  adaptive  since  they  update  the 
solution  using  its  representation  in  a wavelet  basis,  which  concentrates  significant  coefficients  near  singular  behavior. 
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We  used  the  semigroup  method  to  replace  the  nonlinear  partial  differential  equation  with  a nonlinear  integral  equation.  The 
transient  response  of  the  MIS  circuits  has  both  shock-like  as  well  as  smooth  parts.  An  algorithm  based  on  a pseudo-wavelet 
approach,  which  is  automatically  adaptive  to  both  smooth  and  shock-like  behavior  of  the  solution,  was  applied  to  analyze 
this  response.  This  approach  makes  it  possible  to  combine  many  finite-differences,  spectral  and  front-tracking  grid  methods 
into  a collection  of  efficient,  generic  algorithms. 


4.  RESULTS 

Changes  of  incremental  of  hole  density  as  a function  of  time  in  cross  sections  placed  at  different  distances  (x)  from  the 
interface  are  presented  in  Fig.  4.  At  the  interface  and  at  x=0.1pm  increments  are  negative.  At  x=l  pm  the  increments  a the 
beginning  of  the  transient  process  are  negative  and  then,  during  the  equilibriation  period,  are  positive. 


The  inversion  layer  buildap  is  presented  in  Fig.  5,  where  the  move  of  the  zero  current  density  point  as  a function  of  time  is 
depicted.  As  the  gradient  of  the  concentration  of  carriers  in  inversion  layer  rises,  also  rises  the  drift  and  diffusion 
components  of  the  total  current.  This  is  very  crucial  when  calculation  are  based  on  finite  differences.  Wavelets  approach 
enable  to  overcome  these  difficulty. 


Proc.  SPIE  Vol.  4412 


329 


-11  -9  -7  -5  -3  -10  1 lg(t) 

Fig.  5.  Time-dimension  trajectory  of  the  zero-point  current. 

5.  CONCLUSIONS 


The  regions  of  stability  of  the  finite  difference  algorithm  for  Poisson  equation  were  determined.  Approach  based  on 
combined  finite-differences,  spectral  and  front-tracking  grid  methods  into  a cllection  of  efficient,  generic  algorithm  was 
investigated.  Application  of  the  wavelet  approach  increase  the  efficiency  of  the  algorithm,  especially  in  the  equilibration 
period,  when  coupling  of  the  drift  diffusion  equations  through  Poisson  equation  is  very  tight. 
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